Effective Lagrangian for Two-photon and Two-gluon 
Decays of P-wave Heavy Quarkonium Xco,2 smd Xbo,2 states 
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In the traditional non-relativistic bound state calculation, the two-photon decay amplitudes of the 
P-wave Xc0,2 and Xm,2 states depend on the derivative of the wave function at the origin which can 
only be obtained from potential models. However by neglecting the relative quark momenta, the 
decay amplitude can be written as the matrix element of a local heavy quark field operator which 
could be obtained from other processes or computed with QCD sum rules technique or lattice 
simulation. Following the same line as in recent work for the two-photon decays of the S-wave r) c 
and r/b quarkonia, we show that the effective Lagrangian for the two-photon decays of the P-wave 
Xc0,2 and Xbo,2 is given by the heavy quark energy-momentum tensor local operator or its trace, the 
QQ scalar density and that the expression for XcO two-photon and two-gluon decay rate is given 
by the / Xc0 decay constant and is similar to that of r\ a which is given by f v< .. From the existing 
QCD sum rules value for / Xc0 , we get 5 keV for the XcO two-photon width, somewhat larger than 
measurement, but possibly with large uncertainties. 
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I. INTRODUCTION 



With the recent new CLEO measurements [l], H| of 
the two-photon decay rates of the even-parity, P-wave 
++ XcO and 2 ++ Xci states and with renewed inter- 
est in radiative decays of heavy quarkonium states, it 
seems appropriate to have another look at the two- 
photon decay of heavy quarkonium from the standpoint 
of an effective Lagrangian based on local operator ex- 
pansion and heavy-quark spin symmetry, as done for 
the pseudo-scalar heavy quarkonia r\ c and 775 0, H|, for 
which the decay rates for the ground state and excited 
states could be predicted in terms of the J/ ip and T lep- 
tonic widths using Heavy Quark Spin Symmetry (HQSS). 
In the traditional non-relativisitic bound state calcula- 
tion, the two-photon widths for the P-wave quarkonium 
state depend on the derivative at the origin of the spa- 
tial wave function which has to be extracted from po- 
tential models [5j. Though the physics of quarkonium 
decay seems to be better understood within the con- 
ventional framework of QCDQ, unlike the two-photon 
width of <S*-wave r\ c and r]^ quarkonia which can be pre- 
dicted from the corresponding J ftp and T leptonic widths 
using HQSS, there is no similar prediction for the P- 
wave Xc and Xb states and all the existing theoretical 
values for the decay rates are based on potential model 
calculations 0, 0, & ©, M EL M M M EE M E3 • To 
have a prediction for the two-photon width of P-wave 
quarkonia, one need to express the decay amplitude in 
terms of the matrix element of a heavy quark field local 
operator extracted from some known physical processes 
or computed in an essentially model-independent man- 
ner, such as QCD sum rules technique [H, [l9[ or lattice 



simulations [20]. In fact a value of 438 ± 5 ± 6, MeV for 
f Vc and 801 ± 7 ± 5, MeV for f nb consistent with the 
HQSS values of 411 MeV and 836MeV[l i] , respec- 
tively, have been obtained by the lattice group TWQCD 
Collaboration 2l| recently. With similar determinations 
of other quarkonium decay constants, one would be able 
to study QCD radiative corrections and obtain the strong 
a s coupling constant, for example, especially in Xbo,2 two- 
gluon decays where local operator expansion should be a 
better approximation than in Xc0,2 decays. In this paper, 
starting from the two-photon and two-gluon cc — > 77, gg 
and bb — > 77,33 amplitudes, we derive an effective La- 
grangian for the two-photon and two-gluon decays for 
P-wave quarkonium state by neglecting the bound state 
relative quark momenta compared with the large out- 
going photon or gluon momenta. We show that the 
decay amplitude is given by the heavy quark energy- 
momentum tensor which can be obtained from the ma- 



trix element of its trace as < 0|cc|x c q >= 



fxao and 



< 0\bb\xbo >= 



fxbo- We find that the two-photon 



and two-gluon decay rates of Xc0,2 and Xbo,2 are given 
in terms of f Xc0 and f Xb0 , similar to the r/ c and 7/f, two- 
photon decay rates given by / r)c and . 



II. EFFECTIVE LAGRANGIAN FOR X c0,2 

AND xm),2 — > 77 



77 



Following [22, [23] , we consider the amplitude for the 
annihilation of a quark and an antiquark with momentum 
Pi and P2 represented by the diagrams in Fig.fl}: 



A = v{ P2 )(0 1 + 2 )u{ Pl 



(1) 



2 



with 



Oi 



02 = - 



(2) 
(3) 



where (ei,fa) and (£2, fa) are the polarizations and 
momenta of the outgoing photons and mg the heavy 
quark mass. The total energy- momentum of the quark- 
antiquark system is the energy-momentum of the quarko- 
nium bound state defined as Q = p\ + p 2 and mass M. 
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FIG. 1: Diagrams for QQ annihilation to two photons. 



quarkonia with the outgoing photon or gluon having large 
momenta compared to the relative quark-antiquark mo- 
menta as given by the numerator of the amplitude in 
Eq.([7]). We have (rewriting M 2 as 2fa -fa andfa-g£2-fa «i 
as — fa • q^2 ■ fa 4\ in Eq.®): 



M(QQ 77) = ~e 2 Q 



•>,,> A^v(p 2 )T^u(jpi) 
[{fa-fa) 2 /A-m 2 Q ] 2 



c,bl 



(8) 



with A^ u the photon part of the amplitude and the heavy 
quark part T M „ given by: 

A^ v = — 2ej ■ k 2 e 2fi k lL/ + 2t\ ■ e 2 k 2 ^ki L , 

-2e 2 • fa£i M fa„ + (fa ■ fa^ex^v + £2 M £i l /)(9) 
T^ u = (<3V - ?2 M h» (10) 

We see that v(p2)T^ v u(pi) is the matrix element of 

Oqhv = Q(d n - d ^uQ, the heavy quark energy- 
momentum tensor. The photon part can also be written 
in terms of the photon field operator F^, but for simplic- 
ity, we will keep the matrix element form given by A^ v . 
The effective Lagrangian for two-photon and two-gluon 
decay of P-wave Xc0,2 and Xbo,2 states is then given by: 



Using Dirac equation and expanding 0\ and 02 > and 
putting: 

q=pi-p 2 , Q=fa+fa, K = fa-fa (4) 
we have 

n _ 2 n 2 (ei • e 2 (#i- fa) - ie (t2, K, ei, tr)7 <T 7 5 )/2 
Ul ~ CQcb [{pi - kl) 2- m 2 Q } 

_ 2n 2 (-£2 ■ (P2 + fa/2) jx + ei ■ (pi + fa/2) j 2 ) 
G ^ c ' b [(pi-fa) 2 -m 2 Q ] [) 

2n,2 (ei-C2(^2- fa) +it(£2,K,e 1 ,a)j a j 5 )/2 



02 



[(Pi - fa) 2 - TO? 



2^2 (fa ■ (Pi + fa/2) g 7 ! - 61 ■ (p 2 + fa/2) gg) 
[bi-fa) 2 -m 2 ] 



-e 2 g 



(6) 



The P-wave Xc0,2 and Xbo,2 two-photon (two-gluon) de- 
cay amplitudes are given by the P-wave part of the 
QQ ~^ 77? 99 annihilation amplitude which is given by 
the fa ■ q, £± ■ q and £2 ■ q terms in 0i,02. By ne- 
glecting term containing the relative quark momenta q in 
the quark propagator [22| we find, (Q 2 b being the heavy 
quark charge), 

M(QQ^ 77) = ~e 2 Ql b x 



V{P2) 



fa ■ q(—2e\ ■ fa ii + 2ei ■ £2 fa + 2e 2 ■ fa 4\) 



u(pi) 



(fa -fa) 2 /A- m 2 



(7) 



which is now reduced to the matrix element of a local 
operator for two-photon or two-gluon decays of P-wave 



£-ctt{QQ 



77) 
ci 



ic\A llv Q( - d (11) 
e 2 g 2 j6 [(fa-fa) 2 /4-m 2 ? ]- 2 



With the matrix element of Oq^u between the vacuum 
and Xc0,2 or \bo,2 given by (Q 2 — M 2 ): 

< 0\6q^\xo > = T M 2 (- gfiV + QfiQv/M 2 ), 

< 0\e Qllv \ X 2 > = -T 2 M 2 e^. (12) 

where is the polarization tensor for % 2 state, we ob- 
tain the two-photon decay amplitude in a simple manner: 



M(xo 
M(x2 



77) 
77) 



~e 2 Ql 



T A Q 



*[M 2 /4- 



,2 1 12 



T,Ao 



where 



Ao = (-)M 2 (M 2 e 1 



*[M 2 /4- 



£2 - 2ei • fafi 2 • fa) 



m 



2 12 



(13) 
(14) 

(15) 



A 2 = M 2 e llv [M 2 e 1 „e 2v 
+£1 • £2fa M fa^)] 



2(ei • k 2 e 2fl fa„ + e 2 ■ fat ltl k 2v 

(16) 



The above expressions agree with the well-known non- 
relativistic calculation of [5(. The HQSS relation T 2 = 
\/3Tq is obtained by [13] in a calculation of the two- 
photon decays of P-wave quarkonium x,/, J = 0, 2 states 
using the Bcthc-Salpeter wave function and the relativis- 
tic spin projection operators given in this reference and 
in [231 ] which is a precursor of the recent HQSS formu- 
lation of radiative decays of heavy quarkonium [H |3, 
This method allows one to compute the matrix element of 
derivative operator, like the energy-momentum tensor in 



3 



a bound state description of P-wave heavy quarkonium 
state, from which HQSS relations could be obtained. 
However, for QCD sum rules calculation or lattice simula- 
tion, non-derivative operator is simpler to compute. Thus 
instead of working with the energy-momentum Oq^i, op- 
erator, one could work with the trace Bq^ which, by 
applying Dirac equation, becomes a QQ scalar density: 



and 



Qqhh = 2m Q QQ (17) 



v(p2)T^u(pi) = 2m Q v(p 2 )u(p 1 ) (18) 



Then the problem of computing the two-photon or two- 
gluon decays of P-wave quarkonium Xco.2 and Xbo.2 states 
is reduced to computing the decays constants f Xo0 or f Xba 
states, defined as (xo,2 denote here both Xc0,2 and Xbo,2 
states): 



< 0\QQ\XO >= m Xofxo 

Comparing Eq. fHi)) with Eq. (fT9|) , we find 



fxo 



(19) 



(20) 



where we have neglected the binding energy [22j b = 
IrriQ — M and putting vtlq = M/2. This agrees with 
the bound state calculations of [22| and a direct calcu- 
lation of &Qp, v and < 0|QQ|xo > using the expressions 
Eq. (24-26) of [23|. The point we would like to stress here 
is that the local operator expansion allows us to compute 
the two-photon and two-gluon decay amplitudes of Xc0,2 
directly in terms of the f Xa0 decay constant, without us- 
ing the wave function and its derivative at the origin, as 
with that of r) c given in terms of f r]a Q . 

Another quantity of physical interest is the decay con- 
stant f Xl of the P-wave 1 ++ Xci state which enters, for 
example, in B — > Xci^HS Ha] and B — > Xci^H^l de- 
cays. Using expressions Eq. (24-26) in [23[ for xo,i states, 
we find in terms of the derivative of the P-wave spatial 
wave function at the origin 7^(0) : 



■fx a 

fxi 



(87rm Q ) V M 



(87rmQ) 



M 



(21) 



which gives f Xl = 2 fxo- Comparing with the iS-wave 
singlet pseudo-scalar quarkonium decay constant /p[24[ 
(M ~ 2m Q ): 



u = 



32 7rm| 



7^(0) (4mq) 



we have 



f _ 12 ( K(Q) \ f 

/xc0 \K o (0)M J fr >° 



(22) 



(23) 



where TZo{0) is the S'-wave spatial wave function at the 
origin. 

The two-photon decay rates of Xc0.2, Xbo,2 states can 
now be obtained in terms of the decay constant f Xo . We 
find, either by using Eq.pl]) or directly Eq.® for T : 



4iT0 4 a 2 M 3 f 2 



[1 + B (a s /ir)} , (24) 



r77(Xco) ~ (M Xc0 +6) 4 

r 77 (x c2 ) ( m Xc2 + by [1 + B2{asM] (25) 



where B = n 2 /3- 28/9 and Jh = -16/3 are NLO QCD 
radiative corrections [28l I29I l30l|. It is interesting to note 
that, the expression for the Xco two-photon decay rate is 
similar to that for r] c ^: 



r 77 (?7 c ) 



AnQia 2 em M v J 2 r 
(M Vc + b) 2 



a s (20 



(26) 



In the same manner, we have , for the two-gluon decays 
'2\ 47ra 2 M 3 f 2 



T gg(Xco) 



(m Xc0 + by 



■[l + C (a a M], (27) 



r 99 (Xc 2 ) = 77 n 



15 



2 \ 47ra; 



M 3 f 2 



V (M Xc2 +b) 



C 2 (a s /7r)](28) 



where C = 8.77 and €2 = -4.827 are NLO QCD radia- 
tive corrections (2M [29l. |30| . For comparison, the expres- 
sion for Tgg(rj c ) is similar: 



T gg(Vc) = 



^a 2 s M v Jl 
(M Vc + b) 2 



1 



(29) 



We have seen that, the usual expression for the decay 
rate T 77 (xo) m terms of TZ'i{0) is now reduced to the 
simple form Eq. (|24|) by using Eq. (|2Tj) . We note also that 
Eq. (|23|) shows that / Xc0 becomes comparable to f Vc , even 
though, in general, To and f Xc0 are of the order 0(q/M) 
compared with f r)c . 

In the limit of b = 0, the expressions for the two 
decay rates are exactly the same, apart from the de- 
cay constant f xo and f Ve and QCD radiative correction 
terms. The decay rate for Xc2 differs from that of Xco 
only by a HQSS factor. Thus by comparing the expres- 
sion for Xco and r\ c we could already have some estimate 
for the Xco two-photon and two-gluon decay rates. For 
f XcQ of 0(/r) c ), one would expect r 77 (x c o) to be in the 
range of a few keV and that r gg (x c o) is roughly of the 
same size as T gg (r) c ) obtained with the QCD sum rules 
values 18] for the decay constants : f Vc — 374 MeV and 
fxao = 359 MeV which gives, with QCD radiative cor- 
rections (NLO value): r 77 (?7 c ) = 4.33 keV, r 77 (x c o) = 
5.0 keV and r 77 (x C 2) = 0.70 keV to be compared, with 
the latest CLEO result of (2.53 ± 0.37 ± 0.26) keV and 
(0.60 ± 0.06 ± 0.06) keV; and the averages of all cur- 
rent measurements : (2.31 ± 0.10 ± 0.12 keV, 0.51 ± 
0.02 ±0.02) keV, respectively, for Xco and Xc2 two-photon 
width [J. For rj c the prediction from the sum rules value 
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Reference 


r 77 ( X co)(keV) 


r 77 ( X c2)(keV) 


^ _ r 77 (xc2) 
' r 77 (xco) 


Barbieri[5j 


3.5 


0.93 


0.27 


Godfrey [7] 


1.29 


0.46 


0.36 


Barnes [8] 


1.56 


0.56 


0.36 


Bodwin[9] 


6.70 ± 2.80 


0.82 ±0.23 


19+ 015 


GuptaflO] 


6.38 


0.57 


0.09 


Munz[ll] 


1.39 ±0.16 


0.44 ±0.14 


S2+ 016 


Huang[12] 


3.72 ± 1.10 


0.49 ±0.16 


u - lo -0.06 


Ebert[13] 


2.90 


0.50 


0.17 


Schuler[14] 


2.50 


0.28 


0.11 


Crater [15] 


3.34 - 3.96 


0.43 - 0.74 


0.13-0.19 


Wang[16j 


3.78 






Laverty[17] 


1.99 - 2.10 


0.30 - 0.73 


0.14-0.37 


This work 


5.00 


0.70 


0.14 



TABLE I: Potential model predictions for Xc0,2 two-photon 
widths compared with this work. 



experiments. 

For the excited state 2P Xc0,2 states, there has been 
observation of the x' C 2 state above DD threshold, the 
Z(3930) state, at M = (3928 ± 5 ± 2) MeV by_ the 
Belle Collaboration^! which gives T^( x ' c2 ) x B(DD) = 
(0.18±0.05±0.03)keV which implies T 77 (xc 2 ) ^ (0.18± 
0.04) keVH. This would imply r 77 (x' c0 ) ^ (1.30 ± 
0.3) keV, and f x > ~ 195 MeV, comparable with the 
HQSS value of 279 MeV for /^[|. One thus expects 
that T gg (x' c o) in the range 5 - 10 MeV. 

For Xbo,2 the same potential model calculation quoted 
in table I gives the Xf>o,2 two-photon width about 1/10 of 
that for r)b which implies f Xba w (1/3) f Vb , smaller than 
the value obtained from 7£o(0) and 72.^(0) with the Cor- 
nell potential[3l[ which gives f Xb0 = 0A6f Vb . 



III. CONCLUSION 



of f„ c mentioned above is slightly less than the NLO value 
of 5.34 keV obtained with HQSS and is more or less in 
agreement with experiment. Similarly, the prediction for 
XcO from the sum rules value for f Xc0 is however almost 
twice the CLEO value, but possibly with large theoretical 
uncertainties in sum rules calculation for f Xc0 , as to be ex- 
pected. For comparison, we note that the Cornell poten- 
tial model gives f Xc0 — 338MeV[31|. Also a recent QCD 
sum rules calculation^ gives f Xc0 = 510±40 MeV which 
implies a still larger XcO two-photon decay rates. Various 
potential model calculations give r 77 (x c o) hi the range 
1.2 - 6.7 keV and T 77 (x C 2) in the range 0.28 - 0.93 keV 
as shown in table I. From the above expressions for the 
decay rates, the two-photon branching ratios for XcO and 
and rj c would be the same in the absence of QCD radia- 
tive corrections (the two-photon Xc2 branching ratios is 
smaller by 20% with B 1 j/^xs) (X C2 ) = (20 ±1.0)%!). 
With QCD radiative correction, the predicted branch- 
ing ratios are , for a s = 0.28: S 77 (ry c ) = 2.90 x 10~ 4 , 
# 77 (X C0 ) = 3.45 x 10~ 4 , and £ 77 (x c2 ) = 2.55 x 10~ 4 
which are very close to the measured value of (2.4+ - 9 ) x 
10~ 4 , (2.35 ± 0.23) x 10~ 4 , and (2.43 ± 0.18) x 10" 4 , 
respectively^ ■ This shows that QCD radiative correc- 
tions are important in bringing the predictions close to 



By using local operator expansion, we show that the 
two-photon and two-gluon decays of the P-wave heavy 
quarkonium XcO and x&o state can be obtained from the 
heavy quark energy-momentum tensor and its trace, a 
QQ scalar density. The decay rates can then be ex- 
pressed in terms of f Xc0 and f Xb0 decay constant and 
are similar to that of r\ c . Existing sum rules calculation 
for f XcQ however produces a \c0 two-photon width about 
5keV, somewhat bigger than the CLEO measurement, 
but possibly with large theoretical uncertainties. It re- 
mains to be seen whether a better determination of f Xc0 
could bring the Xc0.2 two-photon decay rates closer to 
experiments or higher order QCD radiative corrections 
and large relativistic corrections are needed to explain 
the data. 
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